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Abstract 

We propose a gauge invariant formulation of the effective potential in terms 
of a gauge invariant order parameter, for the Abelian Higgs model. The 
one-loop contribution at zero and finite temperature is computed explicitly, 
and the leading terms in the high temperature expansion are obtained. The 
result is contrasted to the effective potential obtained in several covariant 
gauge-fixing schemes, and the gauge invariant quantities that can be reliably 
extracted from these are identified. It is pointed out that the gauge invariant 
effective potential in the one-loop approximation is complex for all values 
of the order parameter between the maximum and the minimum of the tree 
level potential, both at zero and non-zero temperature. The imaginary part is 
related to long-wavelength instabilities towards phase separation. We study 
the real-time dynamics of initial states in the spinodal region, and relate the 
imaginary part of the effective potential to the growth rate of equal-time 
gauge invariant correlation functions in these states. We conjecture that the 
spinodal instabilities may play a role in non-equilibrium processes inside the 
nucleating bubbles if the transition is first order. 
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I. INTRODUCTION AND MOTIVATION 



In this article we are concerned with the effective potential in gauge theories. It has been 
recognized very early on, that the effective potential is a gauge dependent quantity |l| and 
only a limited amount of information extracted from it is actually physically meaningful. 
This gauge dependence can be understood from several equivalent points of view. The 
effective potential can be identified with the generating functional of 1-particle irreducible 
Green's functions (the effective action) at zero four momentum transfer, and therefore it 
is an off-shell quantity. Alternatively, the effective potential is identified with the energy 
(or free energy) of a particular state (or ensemble) constrained to have an homogeneous 
expectation value of the scalar field. 

The energy, or the free energy, is usually calculated by fixing a particular gauge in the 
path-integral. In a gauge theory, the (complex) scalar fields transform under gauge transfor- 
mations and their expectation value in a gauge fixed state or ensemble is obviously a gauge 
dependent quantity. Despite this shortcoming it has been recognized that certain quantities 
are gauge independent. Dolan and Jackiw |IJ recognized that the critical temperature is a 
gauge invariant quantity and recently Metaxas and Weinberg pi used the Nielsen identities 
@ to prove that the bubble nucleation rate at zero temperature is gauge invariant (we are not 
aware of a similar proof at finite temperature). The gauge invariance of these quantities can 
be understood from the fact that they are associated with homogeneous and inhomogeneous 
extrema of the effective action respectively; these are known to be gauge invariant. 

Considerable effort has been devoted to constructing a gauge invariant effective action 
and effective potential [§J , of which the background field method of Vilkovisky and Dewitt Q 
is the most popular || . These formulations of effective actions are technically formidable and 
do not readily lend themselves to a manageable formulation of equilibrium or non-equilibrium 
descriptions. Furthermore, it has recently been pointed out using the pinch-technique, that 
despite its formal gauge invariance, implementation of the background field method requires 
a gauge fixing parameter for the fluctuations. This leads to a gauge-parameter dependence 
in finite parts of self-energies at finite temperature 0,0], which in turn leads to a gauge 
dependence of the thermal renormalization group beta function as discussed in detail by 
Sasaki |7j. 

Alternative formulations of effective potentials have been offered in terms of a radial and 
angular decomposition of the complex scalar fields |8|,[| or alternatively in terms of gauge 
invariant composite operators [nj. There are several shortcomings in the formulation of 



the effective potential in terms of the radial field variable |8],|9j or composite operators [jTO 
This variable is understood as the "square root" of a composite operator that requires a 
subtraction to be renormalized; shortcomings of this approach had been already recognized 
H . Furthermore in the path integral evaluation there is an ambiguous Jacobian arising from 
the change of variables to radial and angular fields. This Jacobian has to be incorporated 
in the perturbative expansion to obtain a consistently renormalized effective action ||. 

However, even when these technicalities are overcome by some renormalization scheme 
(such as dimensional regularization), it is conceptually unclear how to interpret symmetry 
breaking in terms of the radial field. At the operator level, the radial field variable acquires 
a ground state expectation value even in the symmetric phase as can be seen with a simple 
example of a two-dimensional isotropic harmonic oscillator. Composite bilinear operators 
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typically require subtractions to be renormalized, and their expectation value is therefore 
ambiguous. 

Our motivation is to obtain a gauge invariant description of the effective potential (and 
eventually the effective action) and to use it to provide some preliminary information on 
the dynamics of non-equilibrium processes during the phase transition in gauge theories. 
In view of the above discussion and critique of previous approaches, such an enterprise is 
clearly worthwhile because only a truly gauge invariant description of effective potentials 
can be considered trustworthy in terms of extracting physical quantities such as supercooling 
temperature, latent heat and others that are very important in the quantitative description 
of non-equilibrium features. 

Our program for the construction of an effective potential can be summarized in the 
following steps: a) select the gauge invariant states of the theory, namely, those that are 
annihilated by the first class constraints (such a selection will not involve any gauge fixing), 
b) recognize a gauge invariant order parameter that is invariant under local gauge transfor- 
mations, but transforms non-trivially under the global symmetry that can be spontaneously 
broken, and c) construct the effective potential for this gauge invariant order parameter. An 
attempt to establish a finite temperature framework in terms of gauge invariant states has 
been reported previously [JTTJ] within a different context and with a different goal, but to 
our knowledge it has not been implemented or attempted within the context of the effective 
potential. 

In this article we focus on such a description for the Abelian Higgs model (scalar electro- 
dynamics) and we expect to generalize the procedure and its quantitative implementation 
to Yang-Mills theories in the near future. 

In section II we implement the first step of the program, that is, we select the gauge 
invariant states and order parameter without fixing a gauge, by requiring that the physical 
states be annihilated by the first class constraints of the theory which are recognized as the 
generators of local gauge transformations. Gauge invariant operators are then recognized 
as those that commute with these constraints, out of which we recognize the proper order 
parameter. In section III we explicitly construct the one loop effective potential both at zero 
and non-zero temperature and compare our results with those obtained in popular covariant 
gauges. From this comparison we establish when the gauge-fixed results lead to physical 
(gauge independent) predictions. In this section we also argue that the gauge dependence of 
the usual (gauge fixed) effective potential is not relieved by hard-thermal loop resummation. 
We also provide the high-temperature expansion of the gauge invariant effective potential 
and point out that the "cubic" terms which are typically taken as a signal of the strength 
of a first order transition are in general complex and gauge dependent in fixed gauge path 
integral calculations of the effective potential. 

In section IV we use the gauge invariant effective potential to study the early time 
behavior of spinodal phase separation and the instabilities associated with the spinodal line 
in gauge theories. We establish a correspondence between the imaginary part of the one 
loop gauge invariant effective potential and the rate of growth of correlations in the spinodal 
region. Two appendices are devoted to some technical details. 
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II. THE GAUGE INVARIANT DESCRIPTION: 



The focus of our study is scalar electrodynamics or the Abelian Higgs model whose 
Lagrangian density is 

C = ~F^F„ V + D^D^ - A(0V - (2.1) 

ZV> = {8^ + teA^) (2.2) 

The description in terms of gauge invariant states and operators is best achieved within 
the canonical formulation, which begins with the identification of canonical field variables 
and constraints. These will determine the classical physical phase space and, at the quantum 
level, the physical Hilbert space. 

The canonical momenta conjugate to the scalar and vector fields are given by 

IT = (2.3) 

IT = A i + VM° = -E l (2.4) 

7r + = <fi + ieA°d> (2.5) 

7T = ft - ieA°ft (2.6) 

The Hamiltonian is therefore 

H = J d 3 x j-II • n + tx ] tx + (V0 - ieAcj)) ■ (V0 f + ieAft) + ^(Vx A) 2 + 

A(0 f - tff + A [V • n - ie(vr0 - 7rV f )] } (2.7) 

There are several different manners of quantizing a gauge theory, but the one that exhibits 
the gauge invariant states and operators, originally due to Dirac, begins by recognizing the 
first class constraints (mutually vanishing Poisson brackets). From here there are several 
possibilities: i) the constraints become operators in the quantum theory and are imposed 
onto the physical states, thus defining the physical subspace of the Hilbert space and gauge 
invariant operators, ii) Introduce a gauge, converting the first class system of constraints 
into a second class (with non-zero Poisson brackets between the constraints) and introducing 
Dirac brackets. This is the popular way of dealing with the constraints and leads to the 



usual gauge-fixed path integral representation [12|] in terms of Faddeev- Popov determinants 
and ghosts. 

We will instead proceed with the first possibility that leads to an unambiguous projection 
of the physical states and operators. Such a method has been previously used by James and 



Landshoff within a different context [|T3 ] . 



In Dirac's method of quantization WM there are two first class constraints which are: 



and Gauss' law: 



Q(x,t) = VV -p = (2.9) 
p = ie (dm - ^V) (2.10) 
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with p being the matter (complex scalar) field charge density. 

Gauss' law can be seen to be a constraint in two ways: either because it cannot be 
obtained as a Hamiltonian equation of motion, or because in Dirac's formalism, it is the 
secondary (first class) constraint obtained by requiring that the primary constraint (|2.8|) 
remain constant in time. Quantization is now achieved by imposing the canonical equal- 
time commutation relations 



ir°{x,t),A°(y,t) 

~n\x,t),Ai(y,t) 

7r\x,t),(f)\y,t) 
[7r(x,t),(j)(y,t)] 



-i8(x-y) 

-i5 lj 5(x - y) 

-i5(x-y) 
-iS(x-y) 



(2.11) 

(2.12) 

(2.13) 
(2.14) 



In Dirac's formulation, the projection onto the gauge invariant subspace of the full Hilbert 
space is achieved by imposing the first class constraints onto the states. Physical operators 
are those that commute with the first class constraints. With the above equal-time commu- 
tation relations it is straightforward to see that the unitary operator 



Us 



cxp < i 



n°A + QA 



d x 



(2.15) 



performs the local gauge transformations. Thus the first class constraints are recognized as 
the generators of gauge transformations. In particular, Gauss' law ( |2.9| ) is the generator of 
time independent gauge transformations. Requiring that the physical states be annihilated 
by these constraints is tantamount to selecting the gauge invariant states. Consequently 
operators that commute with the first class constraints are gauge invariant. 

In the Schroedinger representation, in terms of wave-functionals, the canonical momenta 
are represented by hermitian differential operators, and the constraints applied onto the 
states become functional differential equations that the wave-functionals must satisfy: 
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5A (x 



■9[A, 



tfi 
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6A(x) 




[x 



*[A,, 



(2.16) 
(2.17) 



The first equation simply means that the Schroedinger wave-functional does not depend on 
Aq, whereas the second equation means that the wave-functional is only a functional of the 
combination of fields that is annihilated by the Gauss' law functional differential operator. 
It is a simple calculation to prove that the fields 



$(x) 
&(x) 



Ix) exp 



{x) exp 



ic / d 3 yA(y) ■ V y G{y-x) 



-ie / d 3 yA(y) ■ V y G{y-x) 



are annihilated by Gauss' law functional differential equation with G(y — x) 
Green's function that satisfies 



(2.18) 
(2.19) 

the Coulomb 
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V 2 G(y-x)=0 (2.20) 
Furthermore writing the gauge field into transverse and longitudinal components as follows 



A{x) = A L (x) + A T (x) 

V x A L {x) = 

V • A T (x) = 



it is clear that 



5A{x) 



SA L (x) 



(2.21) 
(2.22) 
(2.23) 



(2.24) 



Therefore the "transverse component" At(x) is also annihilated by the Gauss' law operator. 
This analysis shows that the wave- functional solutions of the functional differential equations 
that represent the constraints in the Schroedinger representation are of the form 



(2.25) 



The fields A T ; $ ; & are gauge invariant as they commute with the constraints. The 
canonical momenta conjugate to $ are found to be 



n(x) = n(x) exp 



-ie J d 3 yA{y)-V y G{y-x) 
n f (f) = 7r f (f) exp ie J d 3 yA(y) ■ V y G{y- x) 



^2.26) 
^2.27) 



The momentum canonical to A ; fl is written in terms of "longitudinal" and "transverse" 
components 



n(x) = n,(x) + n T (£) 



(2.28) 



both components are gauge invariant. 

In the physical subspace of gauge invariant wave-functionals, matrix elements of V • II 
can be replaced by matrix elements of the charge density p . Therefore in all matrix elements 
between gauge invariant states (or functionals) one can replace 

U L (x) -> ieV x J d 3 yG(x - y) ($11 - &T^) (y) (2.29) 

Finally in the gauge invariant subspace the Hamiltonian becomes 

H = J d 3 x j-n T • I1 T + n f n + (V$ - ieA T $) ■ (V$ f + ieA T &) + i(V x A T ) 2 + 

A($t$_^ )2 j + ly d 3 y J fapfflGtf-flpfr) (2.30) 

Clearly the Hamiltonian is gauge invariant, and it manifestly has the global U(l) gauge 
symmetry under which $ transforms with a constant phase, II transforms with the opposite 
phase and At is invariant. 
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This Hamiltonian is reminiscent of the Coulomb gauge Hamiltonian, but we emphasize 
that we have not imposed any gauge fixing condition. The formulation is fully gauge invari- 
ant, written in terms of operators that commute with the generators of gauge transformations 
and states that are invariant under these transformations. 

There is a definite advantage in this gauge invariant formulation: the (composite) field 
<&(x) is a candidate for a locally gauge invariant order parameter. The point to stress is 
the following: this operator is invariant under local gauge transformations generated by the 
unitary transformation U\ given by equation (|2.15| ), that is 



U^ffiUX 1 = $ (x) (2.31) 

whereas it transforms as a charged operator under the global gauge transformations generated 
by Q = J d 3 xp(x), that is 

e iaQ $(x)e- iaQ = e iea §{x). (2.32) 

Because the gauge constraints annihilate the physical states and these constraints are 
the generators of local gauge transformations, these states are invariant under the local 
gauge transformations and any operator that is not invariant under these local transforma- 
tions must have zero expectation value. The local gauge symmetry cannot be spontaneously 
broken; this result is widely known in lattice gauge theory as Elitzur's theorem ||15|| . How- 
ever, the global symmetry generated by the charge Q can be spontaneously broken and the 
expectation value of a charged field signals this breakdown. 

From this discussion we clearly see that a trustworthy order parameter must be invariant 
under the local gauge transformations, thus commuting with the gauge constraints, but must 
transform non-trivially under the global gauge transformation generated by the charge. The 
field $ fulfills these criteria and is the natural candidate for an order parameter. 



III. THE EFFECTIVE POTENTIAL 

A. Zero Temperature 

We are now in condition to define the gauge invariant effective potential. Consider the 
gauge invariant state \^;x) such that the expectation value of the gauge invariant order 
parameter <£>(af) in this state is nonzero and space-time constant 

(*;x|$(£)|*;x> = ^ 



(*;xl*;x> 



The effective potential is defined as the minimum of the expectation value of the Hamiltonian 
density in this state, namely 

Veff(x) = n mlD At,.„I,T,.,A ( 3 - 2 ) 



with H being the gauge invariant, Hamiltonian given by equation ( |2.30| ) and Q the spatial 
volume The state |\&; x) is chosen to minimize the expectation value of the Hamiltonian 



subject to the constraint that the expectation value of $ in this state is x- 
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It is convenient to separate the expectation value of $ as 

$(af) = x + v(.%) 



(3.3) 



The one-loop correction (formally of 0(H)) to the effective potential is obtained by keeping 
the quadratic terms in the Hamiltonian. 



H q = n\(\ x \ 2 - 



2\2 



(Vx A) 2 + e 2 A 2 T \ X \ 2 + 



n f n + (VtjXVt^) + 2\rf v (\ x \ 2 - ^ 2 ) + (w f + rfx?} + 

C - J d 3 xd 3 y [U(y) x - nt($ x t] f ) [ n (f) % _ nt(f)x f 



(3.4) 



The transverse components At describe a field with mass m\ = 2e 2 \x\ 2 and only two 
polarizations. The phase of x can be absorbed in II by a global phase transformation under 
which the Hamiltonian is invariant. 

It proves convenient to introduce real fields and canonical momenta as follows 



V 

n 



1 

71 



^(nr-in 2 ) 



(3.5) 
(3.6) 



with {n lj2 ,^i,2} being independent canonical pairs. The non-local part of the Hamiltonian 
is best treated in terms of the Fourier transform of the fields and their canonical momenta, 
in terms of which the quadratic part of the Hamiltonian finally becomes 

H q = m^(|x|) + \ E {^r(k) ■ ft T {-k) + cu 2 (k)A T (k) ■ A T {-k)+ 



^(k^i-k) + uUk^k^-k) + U 2 (k)U 2 (-k) 



ujj.[k) 
k 2 



+ r ]2 (k)r ]2 (-k)u; 2 Jk)\ (3.7) 



where the frequencies are given in terms of the effective masses as 



u T (k) = k 2 + m T 
0Jjj{k) = k 2 + m 2 H 
u; 2 (k) = k 2 + m 2 



m 2 T = 2e 2 \ X \ 2 

m 2 H = 2\(3\x\ 2 -^ 2 ) 

m 2 = 2\(\x\ 2 -fi 2 ). 



(3.8) 
(3.9) 
(3.10) 



The last two terms can be brought to a canonical form by a Bogoliubov transformation. 
Define the new canonical coordinate Q and conjugate momentum P as 



n 2 (fc) 



k 



Ux{k) 



P(k) 



m(k) = ^Q(k) 



(3.11) 
(3.12) 



in terms of which the last term of the Hamiltonian ( |3.7| ) becomes a canonical quadratic form 
with the plasma frequency 
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= = l k2 + 2A(|X| 2 - f)}[k 2 + 2e 2 |x| 2 ]A 2 (3.13) 



There are four physical degrees of freedom. The modes with frequency ujT{k) are the 
two transverse degrees of freedom, the mode with frequency uiuik) is identified with the 
Higgs mode. In absence of electromagnetic interactions (e = 0) the mode with frequency 
ujp(k) represents the Goldstone mode whereas in equilibrium, namely at the minimum of 
the tree level potential, when \x\ = //, it represents the plasma mode which is identified 
as the screened Coulomb interaction, and the transverse and plasma modes all share the 
same mass. However, when the expectation value of the order parameter acquires a non- 
equilibrium value, away from the minimum of the tree level potential (at this order), this 
collective mode does not describe a particle with a Lorentz covariant dispersion relation. 
The frequency clearly shows the combination of the Goldstone dispersion relation and the 
long-range Coulomb interaction typical of a description in terms of the dynamical degrees 
of freedom. 

The lack of manifest Lorentz covariance in the dispersion relation can be understood as 
follows. Although the complex field <fi in the original Lagrangian density is a Lorentz scalar, 
the gauge invariant combination $ given by eq. (|2.18 ) is not, although it is a rotational 



scalar. A particular Lorentz frame has already been chosen in making the transverse and 
longitudinal decomposition of the vector potential ( |2.21| ). The state of lowest available 
energy is expected to be Lorentz invariant, but for arbitrary x there is a constraint in 
the space of functions and these constrained states are not the lowest energy states in the 
functional space. In a scalar theory these states are manifestly Lorentz invariant, simply 
because all fields are Lorentz scalars. With vector fields the situation is more complicated 
and the constrained, gauge invariant states are in general not manifestly Lorentz invariant. 
However the lowest energy equilibrium states (at this order corresponding to x — A*) are 
manifestly Lorentz invariant. We will see in detail in section IV that for x f 1 these states 
are not stationary states of the Hamiltonian, therefore the lack of Lorentz covariance for 
these states is reconciled with their non-equilibrium evolution. 

The quadratic Hamiltonian is now diagonalized in terms of creation and destruction 
operators for the quanta of each harmonic oscillator. The ground state is the vacuum for 
each oscillator and is the state of lowest energy. Therefore the one loop (0(h)) contribution 
to the effective potential is obtained from the zero point energy of the oscillators. Therefore 
accounting for the two polarizations of the transverse components we find: 

Veffdxl) = V d (\x\) + lJ + U H {k) + LO p (k)\ (3.14) 



The normalized wave-functional that satisfies ( |3.1| ) and gives the minimum expectation 
value of the Hamiltonian, thus determining effective potential via (|3.2|) is given by 



V[A T , $ f , $] = iVexp {-- J d 3 x J d 3 yA T {x) ■ A T {y)K T {x - y) 
exp j-^ J d 3 x J (PyrnfflrittifiKHisc - x 
expj-iy d 3 x J d 3 yr]2{x)r] 2 {y)K p (x - y)J (3.15) 
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N = IL 



7P 



K T (x-y) 



K H {x-y) 



d 3 k 
d 3 k 



ik-(x-y) 



ujT(k)e 



uuniky^-^ 



(3.16) 



(3.17) 



(3.18) 



K p (x 



y) 



(3.19) 



(2tt) 3 

d 3 k k 2 U p (k) jk.^-y) 

(2vr) 3 uj.(k) 

This Gaussian wave-functional is clearly gauge invariant, and it has the correct limits: for 
e = (ut = k ;uj p = uo g ) gives the (gauge invariant) wave-functional of free electromagnetism 
times the Gaussian wave-functional of a complex scalar with the U(l) global symmetry 
spontaneously broken, which for \x\ 2 = /i 2 corresponds to the Higgs and a Goldstone mode. 
Writing the fluctuation fields 771 ; r] 2 in terms of $ ; we clearly see that this wave 
functional describes a broken symmetry state since under the global U(l) transformation 
the wavefunctional is changed into an orthogonal wavefunctional in the infinite volume limit. 



The unbroken phase, with x = ; — /r 



m 



> corresponds to the ground state wave- 



functional for free electromagnetism times the ground-state wave-functional of two free real 
scalar fields with equal mass m\pT\. It is straightforward to see that the expectation value 
of the radial variable p = V $t$ is different from zero in this phase and cannot be used as an 
order parameter to signal spontaneous global symmetry breaking as discussed previously. 

The k— integrals in the effective potential (|3.14|) are performed with an upper momentum 
cutoff A. Neglecting a ^-independent term proportional to A 4 as well as terms that vanish 
in the A — > 00 limit, and introducing a renormalization scale k we obtain the following 
(unrenormalized) expression 

A! 
4 



Veff{\x\)=V cl {\ X \) + 



4vr 2 



3m T + m H + m 



1 . 

H In 

16 



4A 2 



3m^ + mjj + m 4 



2m 2 ^m 2 



1 

16 
1 

32 



,fflj,, 

~7T 



,m 



2m^ln(^-)+m^ln(^f) + 



m: 



K 



K 



|) 2 In 



m\, + m 2 + 2Jm 2 m 2 ^ 



K 



3m T + mjj + m 4 , + 6m 2 m 



(3.20) 



The cutoff dependent terms can be absorbed in a renormalization of /i 2 (terms of 0(A 2 ) 
and C(ln(A)) proportional to \x\ 2 ) an d the quartic coupling A (terms of 0(ln(A)) propor- 
tional to |x| 4 )- Using this renormalization prescription we find the following result for the 
renormalized and gauge invariant one-loop effective potential 



eff,R 



(\x\) = ^(\x\ ~ A* ) + ~r~5 \ 00 




3m^ + m H + m A g + 6m 2 m 2 o 



1 

+ 16 



2,4 1„(^) + m < 1„(^£) + (ml - m' T f In 



K 



m. 



m\ + m 2 + 2 



m'rrin 



(3.21) 
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In this expression A and p and all masses are renormalized with the above prescription). 

In the region \\\ 2 < p 2 /3 the Higgs "mass" is purely imaginary , whereas for \\\ 2 < 
p 2 , m? g < 0. Therefore we see that the logarithmic contributions to the effective potential 
from Higgs and plasma modes are imaginary (the last two logarithms in eq. |3.21| ), whereas 
the contribution of the gauge boson is real. The region in which the effective potential 



is imaginary is a region of unstable states fll7|Jl8ll , and the imaginary part of the effective 



potential disguises a non-equilibrium situation whose dynamics will be addressed in section 
IV. This region of instabilities for homogeneous configurations is known as the spinodal 
region. In this region the system is unstable to phase separation, and the imaginary part 
of the effective potential appears as a result of attempting to describe an intrinsically time 
dependent state as a stationary state via analytic continuation. 



B. Finite Temperature: 



The finite temperature effective potential is identified with the free energy density under 
the constraint that the ensemble average of the field is given by a space-time independent 
configuration \. That is 



V T (x) 



X 



—T ln(Tr p) 
Tr $(x)p 
Tr p 



(3.22) 
(3.23) 



with p the ensemble density matrix. In equilibrium and when zero conserved charge is 
considered the density matrix is given by 



p = exp(- — ; 



(3.24) 



In a gauge theory, however, the trace over states in ( |3.22| ) must be defined properly in 
terms of gauge invariant states. Either the physical states are selected and only these are 
used in the trace or alternatively a projection operator must be introduced in the definition 
of the trace |L9| . 



In our approach we select the states as those annihilated by the set of first class con- 
straints, which are therefore gauge invariant as described in the previous section. 

To one-loop order, we have seen that the Hamiltonian is quadratic in terms of gauge 
invariant operators that describe the physical degrees of freedom. Therefore, to this order 
the physical partition function is that of a collection of uncoupled harmonic oscillators for 
each degree of freedom. 

We find the free energy density, which is identified as the finite temperature effective 
potential, to be (/? = l/T) 



F = V eff (\x\;T) = V eff (\ x \;T = o) + 



d 3 k 
pJ (27p 



{2 In 



1-e 



-/3w T (k) 



+ 



In 



-/3w„(k) 



+ In 



(3.25) 



where Kz/dxli ^ = 0) is the zero temperature effective potential given by (|3.14j ), and arises 
from the zero point energy of the oscillators. 
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Just as our gauge invariant approach in terms of gauge invariant operators and func 
tionals allowed us to obtain the ground state constrained wave functional eq. Q3.15|) 



similarly we can obtain the density matrix elements in the Schroedinger representation 
($, & , At\p\A' t , $'). This representation for the density matrix is very useful to study 
non-equilibrium aspects and time dependent phenomena pD| , which is the focus of the next 



section. We find the density matrix elements in the Schroedinger representation to be given 
by 

($, $ f , A T \p\A' T , $' f , = &\pt>\&\ ® (A T \p A \A' T ), (3.26) 

where the density matrices are of the harmonic oscillator type, the full expression is given 
in the appendix I. 



C. Large T expansion: 



The contributions to the free energy from the transverse and Higgs modes are straight- 
forward to obtain by applying the methods developed by Dolan and Jackiw |l| . However the 
contribution from the "plasma" mode is non-standard and requires a more detailed analysis 
which is presented in appendix II. We find the leading high temperature behavior to the 
finite temperature contribution to the effective potential to be given by 



V e ff,T{\x\) 



7T 2 T 4 

~90~ 



rp2 

24 



3m^ + m 2 H + m 2 



T 
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m H m 1 



J^2 



+ 



+ 



(3.27) 



where the dots stand for terms of O(T ) or smaller. Remarkably, the logarithmic terms 
cancel similar terms of the zero temperature part and though this feature is well known in 
the standard cases (with standard dispersion relations for the degrees of freedom), it is a 
new result for the plasma mode. An important feature of this expression is that the terms 
linear in T, which are non-analytic, are complex. Whereas the term from the gauge boson 
mass is real, the terms originating in the Higgs couplings given by the contributions from 
m 2 H , m 2 are purely imaginary on the spinodal regions \x\ 2 < p 2 : /3 for m 2 H and \x\ 2 < p 2 for 
2 These "cubic" terms are usually identified as those responsible for a first order phase 



m 



transition and used to compute quantities relevant to the transition |^T,22,10|. In particular 



these terms determine the supercooling temperature and the latent heat when they are taken 
as the leading indicators for a first order transition. In their study of the electroweak effec- 



tive potential Anderson and Hall [23[ neglected the terms involving the Higgs self-coupling 
keeping only the contributions from the gauge boson and top quark Yukawa couplings which 
are gauge invariant and real to one loop. Arguably such an approximation is justified for 
very weak Higgs couplings. Boyd et.al. P3| recognized that the terms arising from the Higgs 
sector lead to contributions that are imaginary (even after resummation) precisely as pointed 
out above. Thus using these terms to compute the latent heat, supercooling temperatures 
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and even approximate dynamics is at best a crude approximation (even when the imaginary 
parts are ignored) and at worst disguises other non-equilibrium processes which may be 
equally important (see section IV). 



D. Comparison with gauge-fixed results: 



In this section we compare our gauge invariant result to the one loop effective potential 
obtained in the usual standard path-integral representation for several gauge-fixing proce- 



dures. The purpose is to contrast our results with the suggestion of Fukuda and Kugo f25 
that there is a large class of "good gauges" for which the effective potential is gauge in- 
variant. These authors suggested that covariant gauges (including Landau with zero gauge 
parameter), and other specific gauge fixing schemes are such "good" choices. 

In order to make a distinction from the gauge invariant formulation, we write the original 
complex field (j) (not to be confused with the gauge invariant field <&) in the Lagrangian 
density eq. ( |2.1|) as 



{x, t) = -j= [<p R (x, t) + , t) 



(3.28) 



and (p is taken as the (complex) expectation value. To distinguish from the gauge invariant 
case we also introduce the masses 



M 2 g 


= 2A(|yf - 


M\ 


= 2e 2 M 2 


Ml 


= 2A(3M 2 - 



(3.29) 
(3.30) 
(3.31) 



At this point one would be tempted to identify <p with x, because the effective tree level 
masses seem to be the same as in the gauge invariant case under the replacement x ~ * ¥>• 
However we make a distinction between these two expectation values because \ is a truly 
gauge invariant quantity, whereas ip is the expectation value of a gauge transforming field 
in a fixed gauge. 

The one-loop effective potential is given by: 



V = V ^- V f ^ 
1 ^2/3^7 (2tt)3 



In 



2nn 



+ k 2 + M 2 



V 1Q + V 1T 



V 1Q 

v 1T 



i 2 



(0^^ T ^ 

d 3 k 



In 



1 - exp(-pJk 2 + M 2 )) 



(3.32) 
(3.33) 



where the sum are over all particles j with gj degrees of freedom (gj < for ghosts) and 
masses Mj(\<p\), and we have used a result given in [I]. 

The zero temperature contribution is divergent, the k integrals being performed with an 
ultraviolet cutoff A. Discarding a field independent quartic divergence we find the result 
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The finite temperature contribution can be written as 




(3.35) 



The high temperature expansion of I(y) is given by [p]] 

i(y) = Z ^ + T2 y2 ~ l {y2)l ~ ti [Hy2) -\~ c] + ° {y% (3 ' 36) 

where we defined C = 21n(47r) — 27 = 3.9076. The ln(m 2 ) term cancels against the similar 
term in the zero temperature contribution V±o- 



1. Lorentz gauge: 

Dolan and Jackiw ||T| calculate the one loop effective potential both at zero and non-zero 
temperature in Lorentz gauge, which is ghost free, with gauge parameter a. The scalar 
determinant is diagonalized by solving 

ln[A; 4 + M 2 g k 2 + aM 2 M 2 ] = hip 2 + M 2 (a))(k 2 + M 2 (a))} (3.37) 

The resulting masses and effective degrees of freedom are 

1 



Ml{a) 



M 2 ± ^M 9 4 - 4aM|M 2 J ; g ± = 1 (3.38) 

M 9 2 = 2A(M 2 -/x 2 ) ; g g = l (3.39) 

M 2 A = 2e 2 \ V \ 2 ; g A = 3 (3.40) 

Thus the zero temperature part of the one-loop effective potential in Lorentz gauge is 
given by 

1 r d 3 k 

VWM; a) = -J [n H (k) + 3n A (k) + n + (k; a) + n.(k; a)} (3.41) 



J]±(jfe; a ) = Jk 2 + Ml{a) (3.42) 



VLj{k) = y/k 2 + M 2 ;j = H, A. (3.43) 

The gauge dependence (dependence on the gauge parameter a) of the above result is 
explicit. Even for Landau gauge, that is a = 0, expression (|3.41| ) describes five degrees of 
freedom, rather than the four physical degrees of freedom described by the gauge invariant 



result (|3.14j) . With the same renormalization prescription leading to the gauge invariant 



result eq. |3.21| we find in Landau gauge the result 
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V{ip,a = 0) = X(\ip\ 



47T S 



1 

32 



:>>.\r\ + M]j + M, 



+ 



1 

16 



3M 4 ln( !^: ) + M 4 ln( 

K 



Ml 



A M 2 



(3.44) 



which is obviously very different from the gauge invariant result given by eq. |3.21| if (p is 
identified with the gauge invariant order parameter \. 

The effective potential ( |3.41| ) becomes independent of the gauge parameter a for the 
values ip = ; \ip\ 2 = fi 2 . These are the value of the extrema of the tree level potential. 
The gauge dependence appears at 1-loop order and is therefore formally of 0(h) since the 
extrema of the effective potential will acquire 0(h) corrections. We identify the values of 
ip at which the gauge dependence cancels out as the extrema of the effective action to this 
order. Up to an irrelevant constant the gauge invariant effective potential ( |3.14j) and the 
one-loop effective potential in general covariant gauge ( |3.41|) are the same for \x\ 2 = fJ 2 , i.e. 
at the extrema of the effective action. This equality is a consequence of the known result 
that the extrema of the effective action are gauge independent. At zero temperature gauge 
independence at the extrema is also a consequence of the Nielsen identities ||. Therefore we 



find that the suggestion of reference [£5|] is only valid at the extrema of the effective action. 
At finite temperature we find that the 0{T 2 ) is gauge parameter independent and given 

by 



T/(2) 



24 



3M 2 A + M 2 H + M 2 



(3.45) 



This term coincides with that from the gauge invariant effective potential with the identifi- 
cation \ip | = \x\- This contribution determines (to this order) the critical temperature which 
is therefore a gauge independent quantity, whereas the 0(T) and (9(ln(T)) contributions 
are gauge parameter dependent. 

The 0(T) contribution is given by 



T/(l) 



IT 



T 



+ (M 2 H 



+ (A^(a))3 + (M 2 _(a)Y- 



(3.46) 



For a = 0, Ml = M 2 ; M 2 = and this 0(T) term coincides with the 0(T) contribution 
from the gauge invariant effective potential given by eq. ( |3.27| ) if \cp\ is identified with the 
gauge invariant order parameter 

This contribution is of particular importance because it is usually taken as a signal for a 
first order transition and determines its strength, and sometimes used in phenomenological 
equations to describe the dynamics |2l| , |22[| . In the case of first order phase transitions, this 



term is sometimes used to compute the latent heat and the supercooling temperature [10 



Clearly, quantities calculated solely from this term would be physically meaningless because 
of the gauge dependence. Furthermore, this contribution is complex for \(p\ 2 < /i 2 ; only the 
contribution from the gauge boson is real and gauge invariant. Therefore our conclusion 
is that this term is provides the correct gauge invariant 0(T) contribution only in Landau 
(a = 0) gauge and with the identification \<p\ = |%|. However, this equivalence only holds to 
leading order in the high temperature expansion and is not a general feature to all orders, 
as displayed explicitly by higher order finite temperature corrections and also by the zero 
temperature part. 
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2. and gauges: 



In gauge the following gauge fixing and ghost terms are added to the Lagrangian 
density 



1 

"2? 



£fpg = c ] -d 2 - £,e 2 ip(ip + <j> R ) 



(3.47) 



Kastening [j26| describes a useful variant called R% gauge, in which ip and 4>r are treated 
more symmetrically. Its gauge fixing and ghost terms are: 



C 



GF 



1 

'2£ 



FPG 



-«9 2 -£e 2 [(^ + <^) 2 -0 2 ] c (3.48) 



The ghost term is derived as usual by looking at the response of the gauge fixing functional 
under a gauge transformation. The corresponding gauge fixing functional for the gauge 
fixing procedure is f[A^\ = d^A 11 + £e|(</>t<^ ~~ 00)- A piece from the photon determinant 
cancels off half the ghost contribution, so we will treat the ghosts as having g c — — 1. In 
both and R^ gauges we find the following masses and degrees of freedom: 

(3.49) 
(3.50) 
(3.51) 
(3.52) 





M 2 A 


9 a = 3 


Ml 


= £K 


9c = -1 




Ml 


9h — 1 


Mf = M 2 


+ £Ml 


ft = l 



The degree of freedom with M^ is a (gauge-dependent) linear combination of the Gold- 
stone mode and the vector boson. At the tree level minimum \ip\ 2 = [i 2 when M 2 = this 
mode cancels the remaining ghost contribution, leaving a £— independent result, calculable 
just from the photon and the Higgs terms. At the tree level maximum ip = the gauge 
dependence cancels out completely and the ghost decouples. Away from the minimum, 
however, the 1-loop effective potential is £— dependent. 

At finite temperature we find that the 0(T 2 ) term is again gauge parameter independent 
and given by the same expression as in equation ( |3.45| ). 

However the next, 0(T), term is explicitly gauge parameter dependent and given by 



—T 



[3(Ml) 3 / 2 ) - {iM 2 A f 2 + (M 2 H f 2 + (M 2 + £M|) 3 / 2 ] 



(3.53) 



This term depends on £ (unless e = or \ip\ 2 = fi 2 ) and is complex for \cp\ < fi 2 . In Landau 
gauge (£ = 0), it is given by 

— T 



12tt 



[3(M 2 ) 3 / 2 + {M 2 H f 2 + {M 2 f 2 } 



(Landau) 



(3.54) 



and coincides with the 0(T) contribution to the high T expansion of the gauge invariant 
effective potential ( |3.27[ ) if (p is identified with the gauge invariant order parameter \. 

Since to this order only the 0(T 2 ) term enters in the estimate of the critical temperature, 
we see that T c is a gauge invariant quantity (to this order). But comparing the 0(T) term 
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with that in Lorentz gauges should convince the reader that this term is gauge dependent 
and complex in general and one must be very careful in attaching any physical meaning to it, 
such as a criterion for a first order phase transition, its strength and the ensuing latent heat 
and supercooling temperature. Only in Landau gauge do we find that this term is the same 
as in the gauge invariant formulation, although even in this gauge, the higher order finite 
temperature corrections and the zero temperature part are gauge parameter dependent. 



3. Unitary Gauge 



In unitary gauge, even the leading term of Vyr (calculated to 1 loop) is incorrect: ]T 

T 2 



24 



[3(M 2 A )+M 2 H ] 



(3.55) 



In this gauge, higher-loop corrections affect the leading term and must be included in the cal- 
culation to obtain the correct answer (and gauge independent) for the 0(T 2 ) term. Arnold, 
Braaten and Vokos [2~7 showed that a 2-loop calculation restores the correct leading term 



4- Higher-Loop Resummation 

We have seen that the 1-loop effective potential, calculated in various gauge-fixing 
schemes, is explicitly gauge parameter-dependent. One might argue that a resummation 
of higher-loop diagrams would eliminate this dependence. Here we will show that this is not 



the case 28 



Although resummation would naturally be invoked to restore gauge invariance in a partic- 
ular calculation, clearly resummation does nothing to restore gauge parameter-independence 
at T = 0, since higher- loop diagrams are all higher-order in A or e 2 . Thus it is clear that 
the gauge dependence will remain in any calculation that includes the zero temperature 
contribution. 

At finite (but large) T it is a little less obvious, since there is an additional expansion 
parameter: M/T (with M any of the masses). In fact, we have already mentioned that 
resummation in unitary gauge does recover the correct, gauge- independent, T 2 term. We 
now show explicitly in and gauges that the T term remains ^-dependent even after 



resummation of hard thermal loops [28]. 
Expansion Parameters: 

In this section we take e 2 ~ A, and M can represent Mh or Ma- Terms in Vr are 
described as 0(a a f3 hr y c ) with respect to the leading T 2 M 2 term, where [28,29| 



a = AT 2 /m 2 , (3 = XT/m, 7 = y? 2 /T 2 (3.56) 

We will take a » 1 and 7 » 1 but (5 < 1, so that {e 2 ,A} = 0{(3 2 ) and M/T = O(0). 
Our 1-loop high-T expansion is then seen to be an expansion in (3 (not to be confused with 
inverse temperature here). 

We want to resum to 0{(3) (given by daisy diagrams with hard thermal loops). The 
simplest approach is the tadpole method (for a more comprehensive discussion, see ref. 
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), replacing the M/s by thermal masses, which need only be calculated to O(j3 ). We 
then integrate with respect to ip, adding constants to taste, to give Vr- The leading (T 2 ) 
term is unaffected, and remains gauge-independent. 
Thermal Masses: 

The O((3 ) thermal masses are unambiguous and gauge- independent. For the scalars we 
can use the relations Mjj = M\ + V{^((p) and M| = M| + V{ T ((p)/ip (the former is true to 
all orders in gauge). We get 



M 2 H = M 2 H + ^[6e 2 + 8X} : 



Mf = M| + £M 2 + ^[6e 2 + 8A] 



(3.57) 



Only the longitudinal part of the photon gets a thermal (electric screening, plasma) mass 
(given by IIq) to leading order: 



Mr = M\ + 



e 2 T 2 



Ml = M\ 



(3.58) 



With these hard-thermal loop resummed masses, the 0{f3) term in or R^ gauge 
becomes 



-T 



12vr 



2M\ + (1 - £ 3/2 ) M\ + 



g2^2\ 3 / 2 



+ (< + ^[6e 2 + 8A] 



3/2 



+ (M 2 + £Mi + ^[6e 2 + 8A] 



3/2 



(3.59) 



This term is still ^-dependent. We see that resummation (to 0{(3)) does not render the 
effective potential gauge parameter-independent. 



IV. NON-EQUILIBRIUM ASPECTS 

The main motivation for studying an effective potential is to address the issue of sym- 
metry breaking and phase transitions accounting for quantum and thermal corrections. By 
its very definition, the effective potential is an equilibrium quantity because the expectation 
value of the scalar field, that serves as the order parameter, is space-time independent. At 
zero temperature the quantum wave-functional is taken to be a stationary state, at finite 
temperature, the calculation is performed with the equilibrium partition function. Thus 
it is clear that the effective potential is only suitable to describe the equilibrium aspects 
associated with the phase transition. The only equilibrium states correspond to the extrema 
of the effective action which for homogeneous configurations coincide with the extrema of 
the effective potential. Thus quantities like the critical temperature defined as the value 
at which the minima of the effective potential become a maximum, as well as expectation 
values of the scalar field that extremize the effective potential are meaningful quantities that 
are useful to determine whether there is a symmetry breaking phase transition. 

For values of the order parameter away from the extrema, the effective potential is simply 
not a reliable tool to describe the situation and by its very definition it is not meant to be. 
In using the effective potential to address dynamical issues, like "the rolling down" of the 
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order parameter towards the equilibrium configuration the hope is that the time evolution is 
rather slow and the use of the "instantaneous" effective potential is justified as some sort of 
adiabatic approximation. This expectation, however, conceals the physics of the dynamical 
processes during the phase transition 



In second order phase transitions, if the order parameter is initially very small and the 
system is evolving from an initial disordered (high temperature) phase with short range 
correlations towards a final equilibrium broken symmetry state with long range correlations, 
the dynamical process is that of phase separation and growth of correlated domains. This 
dynamics is not captured solely by the time evolution of the expectation value of the scalar 
field, but the growth of fluctuations will be manifest in the time dependence of the correlation 
functions of this field. During the early stages of a second order phase transition, when the 
order parameter is very small (near the maximum of the effective potential) small amplitude, 
long wavelength fluctuations become unstable and grow []30[| . The two-point correlation 
functions of the scalar field reflect this instability and grow exponentially, as the order 
parameter "rolls down" the potential hill. 

If the phase transition is of first order, then there are free energy barriers that the 
system has to overcome to reach the equilibrium state and the phase transition is driven by 
nucleation, in which large amplitude configurations become unstable and grow. However, if 
the first order transition is very weak (that is to say that there is a substantial amount of 
phase mixing) nucleation and phase separation occur on similar spatial and time scales and 
the long-wavelength instabilities will still be important. 

The information on thermodynamic instabilities that lead to phase separation is con- 
tained in the effective potential both at zero and finite temperature in the form of an imag- 
inary part. 

Weinberg and Wu |T7[ have shown in a beautiful paper, that the imaginary part of 
the effective potential in scalar theories determines the growth rate of the scalar field two- 
point correlation function for Gaussian states centered at zero expectation value (in field 
space). Boyanovsky and de Vega |Tj| have studied how this growth of correlations affects the 
time evolution of the order parameter and concluded that the use of the effective potential to 
describe the time evolution of the expectation value of the scalar field is not only unwarranted 
but completely misleading and unreliable whenever the initial value of the expectation value 
is in the "spinodal region" (the region in the tree-level potential where the second derivative 
is negative). 



A. The Imaginary Part 

We notice that in the "classical spinodal" region, where m 2 H < or < \x\ < n/VS, 
the Higgs mode has a band of unstable wave-vectors with imaginary frequencies for k 2 < 
fi 2 — 3\x\ 2 - The plasma mode has a band of unstable wave vectors with imaginary frequencies 
for k 2 < fi 2 — \x\ 2 in the spinodal region < |x| 2 < /i 2 - This new spinodal line ranges from 
the maximum to the minima of the tree level potential. 

When the gauge coupling is switched off, in the absence of long-range forces, this is 
recognized as the "spinodal" region for the would be Goldstone modes PU[. However in 



the presence of long range forces, the instabilities are much more severe for long- wavelength 
fluctuations, as can be seen from the infrared behavior of the plasma frequency eq. ( 3.13| ). 
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As a consequence of these unstable modes, both the zero and finite temperature effec- 
tive potential acquire an imaginary part for all values of \x\ between the maximum and 
the minima of the tree level potential. Therefore we see that unlike the case of a scalar 
order parameter in which there is a thermodynamically stable region in the phase diagram 
(between the classical spinodal and the coexistence line), in this case the unstable region 
covers the whole of the phase diagram below the coexistence curve. 

Thus we here obtain one of the important conclusions of this work: the gauge invariant 
effective potential to one-loop order is complex from the maximum to the minimum of 
the tree-level potential. There are no homogeneous configurations corresponding to quasi- 
equilibrium states away from the minimum. Explicitly, the imaginary part of the gauge 
invariant effective potential at zero temperature is given by: 

rAl(lxl) 




Im(V eff )= ±4^2 y k z dk^/kl(\x\)-k* 9(^-3|x| 

+ J* * kdk^kl{\x\) - k 2 ^2e 2 \ X \ 2 + k 2 Q{^ 2 - 

kldxl) = 2A(^ 2 - 3| X | 2 ); *2(|x|) = 2A(^ 2 - | X | 2 ) (4.1) 

where the ± is determined by the direction of the analytic continuation in the frequencies. 

The presence of an imaginary part of the effective potential has been sometimes justified 
as a failure of the loop expansion and "corrected" by the Maxwell construction which yields 
a convex free energy. However, just as the Van der Walls loop in the equation of state 
for liquid-gas systems is a signal of thermodynamic instabilities, the imaginary part of the 
effective potential signals the presence of similar instabilities in the quantum system. 

The flat region of the Maxwell constructed free energy indicates that the system will be 
found in a coexistence of phases but offers no information on the non-equilibrium processes 
leading to phase separation. In this respect the effective potential with its imaginary part at 
least provides some restricted information on dynamical processes. As discussed in references 
|T71 . |T8H this imaginary part determines the decay rate of unstable Gaussian states. 

In the one-loop approximation, the Hamiltonian for the modes with wavevectors in the 
unstable bands correspond to inverted harmonic oscillators for which the analysis of ref- 
|IT| , |i^ , |2Ti| , |3"Tf can be applied. The effective potential ( |3.14| ) is complex because the 



erences 



modes corresponding to inverted harmonic oscillators were treated as ordinary harmonic 
oscillators, i.e. by analytically continuing the zero point energy for these oscillators. The 
Gaussian wavefunctionals and density matrix have complex kernels for the fields 77^2 reflect- 
ing this analytic continuation. 

In terms of the shifted fields 771 , 772 with zero expectation value, the Gaussian wave- 
functional (|3.15| ) and the density matrix (see appendix I) are "centered" at the origin in 



field space, and for the wavevectors in the unstable bands, they are not stationary states of 
the Hamiltonian. 

In this one-loop approximation we can use the results of previous investigations of similar 
situations |T7|Jl~8|j20|j3T|| to study the early time evolution of initially prepared non-equilibrium 



states. 
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B. Early time evolution of non-equilibrium states 



From the previous discussion it is clear that the imaginary part of the effective potential 
conceals a dynamical non-equilibrium situation associated with the instabilities towards 
phase separation. This is a dynamical situation that must be treated as a time dependent 
problem. The quantum states that lead to the effective potential (and free energy) are not 
stationary states. 

Since only the Higgs and plasma modes have unstable bands we will focus only on these. 
Consider the situation in which an initial Gaussian density matrix has been prepared such 
that the gauge invariant scalar field $ has initial expectation value in this ensemble given 
by Xi an d that the kernels have positive frequencies. That is, the density matrix in the 
Schroedinger representation is given by the expression in appendix I but with the following 
real kernels: 



K ( H ] {k,t = 0) = ~W H {k) coth 



1 



K^(k,t = 0) = --Wp(k) coth 
Wjj(k) > , Wp{k) > 



W H (k) 
T 

W P {k) 
T 



K§\k,t = 0) 

K%\k,t = 0): 



W H (k) 



W H (k) 
T 



sinh 
W P (k) 



sinh 



W P (k) 
T 



(4.2) 

(4.3) 
(4.4) 



The frequencies Wn{k) , Wp{k) determine the initial conditions on the state as discussed 
in refs. |l7yi8i^ |3l] . We can now use the results of references [|T7Hl8tp0i^ll to determine 
the time evolution of this initial state. The early time evolution of the expectation value 
X will be determined by the classical equations of motion with small quantum and thermal 
corrections. However at longer times the growth of the unstable modes makes the one-loop 
approximation to the dynamics unreliable as in the problem of domain growth in scalar 
theories |2(]|JIJ]. Clearly, the most important dynamics is described by the evolution of 
the unstable modes and determined by the time evolution of the kernels. We now use the 
method discussed in reference EDI] to obtain the following time evolution of the kernels: 



K$\k,t) 
Kg\k,t) 



K^\k,t = 0) d\n(\U k (t)\ 2 



dt 



K§\k,t = 0) 
\U k (t)\ 2 

KU(k,t = 0) . dM\V k (t)f 

\v k (t)\ 2 1 

K( 2 \k,t = 0) 



dt 



(4.5) 
(4.6) 
(4.7) 
(4.8) 



IW)| 2 

where the mode functions £4(t) = Uik(t) + iU2k(t) , Vik(t) + iV 2 k(t) obey the following 



evolution equations |2C 
d 2 



dt 



2 +k 2 + 2\(3\ Xc i\ 2 (t)-ix 2 ) 



U k (t) = 
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U lk (0) = 1 , Z7 2Jfe (0) = , U lk (0) = , U 2k (0) = W H (k) 

d 2 



dt 2 + 



k 2 + 2X(\ Xd \ 2 (t) - ^ 2 )j [fc 2 + 2e 2 \ Xcl \ 2 (t)\ jk 2 
V lk {0) = 1 , ^(0) = o , K lfe (0) = , F 2fc (0) = W P (k) 



V k (t) = 



(4.9) 



(4.10) 



where Xd is the classical evolution of the "zero mode" |p0| . Obviously the mode functions 
U k {t) , V k {t) grow almost exponentially for wave- vectors in the respective unstable bands 
when is m the corresponding spinodal regions. 

Consider the equal-time two-point correlation functions of the gauge invariant operator 



V, V 1 



d 3 k 



D ik-(x-y) 



\u k {t)\ 2 

2W H (k) 



coth 



W H (k) 



+ 



2W P (k) 



coth 



W P (k) 



(4.11) 



When Xd(t) remains in the spinodal, it is clear that the growth rate of (77(0, t)rf(0, t)) is 
related to the imaginary part of the effective potential given by eq. ( [4. 1| ) at zero temperature. 

At this point we recognize an important payoff of the gauge invariant description. In 
order to compute gauge invariant correlation functions from the gauge variant operator 
4>(x, t) we would have to append a line integral of the (time dependent) gauge field, with the 
ensuing path ambiguities and complications. The formulation in terms of gauge invariant 
order parameters from the start overcomes these difficulties and allows to extract physically 
meaningful correlation functions that provide dynamical information on non-equilibrium 
processes. 

Just as in the case of a scalar field theory undergoing spinodal decomposition 
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equal-time correlation function will grow exponentially at early times because of the unstable 



modes. This growth of is the hallmark of the dynamics of phase separation [17,18,30 



If the phase transition is weakly first order in the sense that a considerable amount 
of phase separation and mixing occurs during nucleation, then the growth of correlations 
described by the dynamics of the unstable modes will be also a dominant process. 

As noted before, unlike the case of a scalar field transforming under a discrete symmetry, 
in this case the spinodal region reaches all the way to the minimum of the potential (at 
least in the one-loop approximation) and the long- wavelength instabilities are enhanced by 
the long-range forces. This feature provides interesting possibilities if the phase transition 
is strongly first order and nucleation is the dominant mechanism. Typically after nucleation 
of a critical bubble, the order parameter inside the bubble is very close to its equilibrium 
value. If this value is smaller in magnitude to that of the equilibrium configuration, spinodal 
instabilities may be still present inside the bubble. 

If the critical droplet has rather small radius it will have to grow to a size ~ l/(|x| — /-0? 
with \x\ the value of the order parameter inside the bubble, before any instability develops 
inside the bubble. Therefore the fact that the spinodal line reaches to the minimum may lead 
to the tantalizing possibility that for large critical bubbles there may still be non-equilibrium 
processes inside the bubble if the order parameter inside is smaller in magnitude than the 
equilibrium value. 

Clearly all these possibilities will have to be studied in deeper detail and we expect a 
strong dependence on the values of the gauge and Yukawa couplings (in the case of fermions). 



22 



Furthermore our analysis of the time evolution of the spinodal instabilities only holds at very 
early times after the preparation of the initial state. At longer times a non-perturbative 
gauge invariant scheme will have to be used to determine the dynamics. We hope to be able 
to implement the variational schemes of reference to the gauge invariant formulation. 



V. CONCLUSIONS 

In this article we have presented a formulation of the effective potential in terms of a 
gauge invariant order parameter for the case of the Abelian Higgs model. The gauge invariant 
states of the theory are those annihilated by the first class constraints and gauge invariant 
operators are those that commute with these constraints. We recognized an order parameter 
that is invariant under the local gauge transformations but transforms as a charged operator 
under global phase rotations, its expectation value in the lowest energy state therefore signals 
the breakdown of the global U(l) symmetry. 

We evaluated the one-loop contribution to this gauge invariant effective potential both 
at zero and non-zero temperature and obtained its high temperature expansion. We found 
that the effective potential is complex both at zero and non-zero temperature and that the 
spinodal (thermodynamically unstable) region extends from the maximum to the minimum 
of the tree level potential (in this approximation). The gauge invariant effective potential 
was compared to the effective potential obtained in several covariant gauges, we found that 
the dependence on gauge parameter cancels only at the extrema, but for all other values 
of the order parameter the gauge fixed effective potentials are gauge parameter dependent 
at zero and non-zero temperature. In particular in a high temperature expansion only the 
0(T 2 ) contribution is gauge invariant, whereas the 0(T) , (9(ln(T)) ■ • ■ depend on the gauge 
parameter and are complex. In general these cannot be taken as trustworthy quantities to 
extract information on the strength of the transition and its features, such as latent heat 
and supercooling temperature, although the contribution of the gauge boson masses (we did 
not study fermions) are gauge independent and real to this order. 

The imaginary part of the one-loop effective potential determines the spinodal lines, 
which identify the region of long- wavelength instabilities. Unlike the case of a scalar field 
with discrete symmetries, the spinodal line encompasses all values of the homogeneous order 
parameter ranging between the maximum and the minimum of the tree level potential in 
the one-loop approximation. 

Non-equilibrium dynamical aspects were then studied in terms of gauge invariant correla- 
tion functions. It is in these dynamical situations out of equilibrium that the gauge invariant 
formulation has the greatest impact, since in general, correlation functions involve non-local 
line integrals of the gauge field. We obtained the early time behavior of the gauge invariant 
correlation functions for initial conditions in which the order parameter is in the spinodal 
region, these correlation functions grow exponentially at early times as a consequence of the 
instabilities. 

We conjectured that these spinodal instabilities may play an important role in first order 
phase transitions, in that they may be responsible for non-equilibrium dynamics inside the 
nucleating bubbles. Clearly this possibility will have to be studied further. 

There are several important avenues to pursue: higher order calculations and the imple- 
mentation of variational or Hartree-like approximations in the gauge invariant formulation to 
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address the long-time behavior and to obtain a more solid understanding of non-equilibrium 
processes. Furthermore, it is important to generalize the resummation program to the gauge 
invariant effective potential, and to generalize these new methods and results to non-Abelian 
gauge theories in particular the electroweak theory. Study on these issues is underway. 



VI. APPENDIX I 



As shown in section II, in the one-loop approximation the Hamiltonian becomes a sum of 
independent harmonic oscillators for the gauge invariant normal modes and in this case the 
explicit form of the density matrix in the Schroedinger representation is available (see 



and references therein). This form of the density matrix is particularly convenient to study 
real-time dynamics, since the Liouville equation becomes a functional differential equation, 
which in this case can be solved explicitly for the time dependence [^C] . 

The equilibrium density matrix that describes the system at temperature T — is 
given in the Schroedinger representation by: 



($, A T \p\A' T , $ /f , = N[T] x 



exp^ ld 3 x \d 3 y A T (x) ■ A T (y) + A' T (x) ■ A' T (y) K$\x - y) + A' T (x) ■ A T (y)K ( T Z) (x - y) 



-(2), 



x exp \ \d 3 x ld 3 y [r]i(x)r] 1 (y) + r]'i(xWi 



K^'ix-y) +r][(x)7] 1 (y)K ( H > (x - y) 



(2), 



x exp( fd 3 x [d 3 y [r] 2 (x)r] 2 (y) + rj' 2 (x)r] 2 (y)} (x - y) + i 2 (x)r] 2 (y^Kf\x - y)\ (6.1) 
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Tu 2 T (k) 

_ ik-(x-y) 
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(2tt) 3 



u^(k) sinh 



k 2 uj p (k) 
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(6.2) 
(6.3) 
(6.4) 

(6.5) 
(6.6) 

(6.7) 
(6.8) 



The trace of this density matrix gives the gauge invariant one-loop partition function, 
whose logarithm (divided by —f3) gives the one-loop effective potential at finite temperature. 
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VII. APPENDIX II 



In this appendix we provide the essential ingredients to obtain the high temperature ex- 
pansion of the contribution of the plasma mode to the finite temperature effective potential. 

The first two terms in the free energy density given by eq. ( |3.25| ) (the transverse and 
Higgs modes) are of the usual form with dispersion relations of the form u(k) = y/k 2 + m 2 
and their asymptotic high temperature expansion can be obtained by following the steps 
described in Jj]. The third term requires special attention because of the unusual dispersion 
relation of the plasma mode. This term can be written as 



d 3 k 



111 



-f3ui p (k) 



rp4 



x dxln 



o 

m 



Vl — mo 



£ 
J^2 



vl 



J^2 



The function I(yf,y 2 ) can be written as 

d d 



<](,dh fafb I{ - a ' h) + /(0 ' ^ + ^ 0) " /(0 ' 0) 



(7.1) 



(7.2) 



io Jo 

The last three terms are of standard form and can be calculated following the procedures of 
reference M. Using the identity (3.11) of reference [IJ we find 

rA 



d d T , 1 d d 

I (a, b) = 

da db 2 da db 



J xdx\J (x 2 + a)(x 2 + b) + 



rA 00 
/ x 2 dxy^ 



2nlx 



(x 2 + a)(x 2 + b) + (2nlx) 



(7.3) 



where we introduced an upper momentum cutoff to regulate the individual terms. The 
integral in the first term is the same as the plasma contribution to the zero temperature 
effective potential in eq. ( 3.14Q but in terms of the variables a, b. This integral has a 
logarithmic cutoff dependence. The second term does not have infrared divergences for 
a, b = and can be expanded in power series of a, b. In this series, the first term, with 
a = b = has a logarithmic cutoff dependence that exactly cancels (upon integration 
on a, b) that of the first term. The remaining series is both ultraviolet and infrared safe. 
Finally only the finite terms from the zero temperature effective potential but in terms of the 
variables y\\ y\ contribute to I(yf, y%) and we obtain the result quoted in equation (|3.27|). 
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